Surface spin-flop phases and bulk discommensurations in antiferromagnets 
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Phase diagrams as a function of anisotropy D and magnetic field // are obtained for discommen- 
surations and surface states for a model antiferromagnet in which H is parallel to the easy axis. 
The surface spin-flop phase exists for all D. We show that there is a region where the penetration 
length of the surface spin-flop phase diverges. Introducing a discommensuration of even length then 
becomes preferable to reconstructing the surface. The results are used to clarify and correct previous 
studies in which discommensurations have been confused with genuine surface spin-flop states. 



Recent experimental results on Fe/Cr(211) superlattices jl],g] have stimulated theoretical and numerical work 
aimed at explaining the occurrence of series of phase transitions in samples with an even number of Fe blocks. 
The antiferromagnetically-coupled Fe blocks can be described phenomenologically by a chain of classical X-Y spins 
with antiferromagnetic interactions and uniaxial spin anisotropy D. For a chain of infinite length, application of an 
external magnetic field H in a direction parallel to the easy axis (the spin direction in zero field) results in a transi- 
tion at a finite field H = Hsf(D) to a spin-flop phase in which the spins on the two sublattices are approximately 
perpendicular to the applied field, and nearly opposite to each other g|. 

For finite chains the situation is more complex. In 1968 Mills Q suggested that spins in a layer near a free surface of 
such an antiferromagnet could rotate into the flopped state at a field H' SF significantly less than H$f- To understand 
this further we present numerical and analytic studies for an antiferromagnetic chain for all values of H and D. 
We show that a consideration of both surface states and discommensurations is necessary to fully understand the 
phase diagram. It transpires that discommensurations, while they are not surface phenomena per se can provide the 
minimum energy configuration for both semi-infinite and finite systems. The results allow us to unify certain aspects 
of previous work |l],|2|,|]-[7| while correcting others. 

Consider a chain of classical spins described by a Hamiltonian 



n = J2 



Oi + \) — HcosOi — ^cos(2#i) 



(1) 



where the exchange coefficient has been taken as the unit of energy, 9i is the angle between the zth spin and the 
direction of the field H, and D is a two- fold spin anisotropy. One can think of Qi as the direction of the magnetization 
in the ith layer of an antiferromagnet, or in the ith Fe layer in a Fe/Cr superlattice. It is appropriate to work at zero 
temperature in one dimension because in-plane fluctuations are not important for the superlattice physics. 

For a bulk system, in which the sum in (Q) runs from — oo to +oo, the ground state is either ferromagnetic (F) 
with all spins parallel to the field, antiferromagnetic (AF) with the spins alternating between and tt, parallel and 
antiparallel to the field, or a spin flop (SF) phase in which the spins alternate between +(f> and —6, for some <fi in 
the range < < n/2. The regions where these different phases are stable are indicated in Fig. H, where the AF 
phase is above the line D(A — D) = H 2 and to the left of H = 2, while the F and SF phases are separated by the line 
D + H = 4 The additional structure in the AF region refers to discommensurations and surface phases. 

We first describe the phase diagram of a chain which is constrained by suitable boundary conditions to include a 
discommensuration ||. Note that the bulk AF ground state is degenerate with two possibilities: 0* equal to for i 
even and tt for i odd, or vice versa. A discommensuration is a boundary or interface between two regions corresponding 
to these two possibilities; in particular, a configuration in which Qi n tends to and 02n+i to tt as n tends to +oo, 
while 6*2n tends to tt and #2n+i to as n tends to — oo. For all values of D above the dashed line in Fig. [|, the 
discommensuration of minimum energy is of the Ising type, . . . 0, tt, 0, tt, 0, 0, 7T, 0, it, . . . , with two adjacent spins 
parallel to each other and to the field, in the middle of what is otherwise an antiferromagnetic configuration. 

In the AF region below the dashed line in Fig. [j], "flopped" discommensurations of different length have lower 
energies than the Ising discommensuration. A flopped discommensuration of type (2m) consists of a "core" of 2m 
spins in which the spin configuration resembles that in a bulk SF phase, located between "tails" , each of which rapidly 
reverts to the configuration of the corresponding AF phase with increasing distance from the core. See, for example, 
Fig. 1 in jrofl . As D decreases, the distinction between the "tails" and the "core" is less clear, but we continue to 
use the same label for the discommensuration which evolves continuously from (2m) at larger D. One can think of 
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such a discommensuration as composed of a pair of AF-SF and SF-AF interfaces centered at the points where the 
corresponding tail joins the core. At low values of H, the discommensuration (2) has the lowest energy, but upon 
approaching the bulk AF:SF phase boundary, one finds — see Fig. || for details omitted from Fig. [l] — a sequence of 
phase transitions to (4), (6), . . . as H increases. Our numerical procedures, which used effective potential techniques 
p| |l3| to construct surface states by extrapolating from the bulk E3], found values of 2m up to 14. We were able to 
trace the first-order lines separating the different (2m) phases down to a value of D between 0.1 and 0.4. We found 
no evidence to suggest that these lines end in critical points, and we believe it likely that they will persist all the way 
down to D = 0. 

These transitions reflect the discrete nature of the spin chain. Thus it is not surprising that they are absent in 
the continuum approximation employed in fl(i| for small values of D. That study showed that the width of the 
discommensuration tends to oc at the AF:SF boundary, in agreement with what we found for larger values of D. 

The boundary between (2) and the Ising discommensuration is second order, and represents the limit of stability of 
the latter as D decreases. An analytic calculation yields the equation (D + H — 1) _1 = 5/3 + D — H, in good agreement 
with our numerical calculations, and those in |hJ when H = 0. The boundaries between the Ising discommensuration 
and (4), (6) etc., are first order, and were obtained numerically. The triple points at which the phases AF2, (2m) 
and (2m + 2) meet tend to an accumulation point, Q, located at H w 1.58,13 ~ 0.78. Presumably this is where 
the energy to create a pair of AF-SF and SF-AF interfaces (infinitely far apart) is equal to the energy of an Ising 
discommensuration. 

We now consider the surface states of a semi-infinite chain. If a ground-state configuration of an infinite chain is cut 
in two while the spins are held fixed, there results two ideal or unreconstructed surfaces. Sometimes the energy of the 
surface state can decrease through a local reconstruction in which spins near the surface are altered from their bulk 
values by amounts which tend asymptotically to zero with increasing distance from the surface. The energy change 
during such a reconstruction can be calculated from the appropriate Hamiltonian, which is ([!]) with the sum going 
from i = to +00, even though the total energy is not well defined. Hence the reconstructed surface of minimum 
energy is (usually) well-defined. 

While the surface reconstruction in the spin-flop phase occurs smoothly, so that there are no phase transitions, the 
situation in the case of antiferromagnetic surface states is more complicated. There are two types of surface states, 
A and B, with unreconstructed versions having the surface spin parallel (#0 = 0) or opposite (0q = tt) to the field 
direction respectively. Even after reconstruction takes place, so that 6q has changed, the A-type (-B-type) surface can 
still be identified through the fact that (9 2ra tends to (tt) as n — > 00. 

Throughout the AF region, A-type surfaces do not reconstruct. The behavior of B-type surfaces is more complicated. 
In region AFi in Fig. [|, the unreconstructed surface has the lowest energy. In region AF2, which meets AFi along a 
line H = 1 for arbitrarily large D, there is a set of degenerate (equal minimum energy) reconstructed surfaces, 



where \2n) consists of 2n spins in an antiferromagnetic arrangement, followed by two spins parallel to the field, and 
then the bulk antiferromagnetic phase. One can think of this reconstructed surface as an Ising discommensuration 
located a certain distance from the surface. Because the "tails" of this discommensuration have zero length, it does 
not interact with the surface, and its energy is independent of its distance from the surface as long as 9q — 0. (This 
artificial degeneracy would be lifted if weak longer-range interactions were included in (jl]).) 

In the AF3 region in Fig. [|, the B-type surface again reconstructs, but the spins near the surface are no longer locked 
in Ising positions. Just as in the AF2 region, one can think of the surface state as consisting of a discommensuration 
located a finite distance from the surface, but now this discommensuration is of the flopped type with a core of length 
two, and tails extending out on either side of the core. Because of the tails, the discommensuration interacts with 
the surface, and the minimum surface energy occurs when it is some specific distance 2n from the surface, depending 
on D and H. Thus in AF3 one finds genuine spin-flop surface states. As H increases, the discommensuration moves 
further from the surface. It does this, at least when D is large, discontinuously in steps of 2, via a series of first-order 
phase transitions, some of which are shown in Fig. ||, where they extend leftwards from the point P. Numerically we 
have seen states with 2n up to 14, and our results are consistent with n tending to infinity at the right side of the 
AF3 region, which our analytic calculations, in agreement with ||, show to be the line 



The upper boundary of the AF3 region extending from O to P is a continuous (second-order) transition representing 
the limit of stability of the Ising surface phase [0), see (|2]), as D decreases. An analytic calculation yields the implicit 
equation 



[0) = 007r07r07T..., 



[2) = OttOOttOtt. . . , 



(2) 



D = y/l + H 2 - 1. 



(3) 



(2 + D-H-l/a) 
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2 + D + H -a, 



(4) 
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where a = H + D + 1/(1 — H — D). Thus the point P. where all the phases \2n) come together, lies at H = 4/3, 
D = 2/3, the intersection of @ and (@). Both (|) and (§) agree with our numerical results. 

We find that the first-order lines extending downwards and leftwards from P in Fig. 0, separating the canted phases 
[2n) from [2n + 2), end in critical points as D decreases. The larger the value of n, the further the line extends 
towards the origin, but presumably for any finite value of n the difference between the phases [2n) and [2n + 2) 
eventually disappears at some finite value of D. Because this value decreases with increasing n, it is plausible that 
the corresponding critical points accumulate at the origin. 

As is evident in Fig. |], the region AF3 becomes extremely narrow as D decreases. The left boundary approaches a 
parabola D = 0.5H 2 to within numerical precision, which is asymptotically the same as (||). We nonetheless believe 
that the width of AF3 remains finite as long as D > 0. Numerical evidence for this is that the value of the surface 
spin 6*o, at the left edge of the AF 3 region (that is for H just large enough to produce the surface spin- flop phase) 
tends to a value near 7r/3 as D goes to zero. This indicates that even for very small D the discommensuration at the 
threshold field is still a finite distance from the surface. 

Between AF3 and the AF:SF bulk phase boundary lies region AF4 in which the flopped discommensuration is 
repelled by the surface, so that its minimum energy location is in the bulk infinitely far away from the surface, as 
noted in ref. ||. Thus there is not a minimum-energy reconstructed B surface, or, properly speaking, a "surface 
spin- flop phase" in region AF4. It would probably be better to identify AF4, thought of as part of the _B-type surface 
phase diagram, as the "discommensuration phase" , since the minimum energy surface will always be of the A-type. 

In retrospect it seems likely that the broadening of the SSF transition mentioned in the abstract of ref. || actually 
refers to broadening of the bulk discommensuration which, as noted above, occurs as H approaches the AF:SF 
boundary inside region AF 4 . It seems that no earlier work has correctly identified the stable SSF phase at small 
values of D, characterized when it first appears with increasing H by a surface spin with a value very near 7r/3. The 
narrowness of the AF3 region for small D may be why it was overlooked. 

For a finite chain of length L, the minimum energy is achieved with A-type surfaces at both ends when L is odd. 
This is not possible if L is even, so when H is small there will be an A-type surface at one end and a J5-type at the 
other. As H increases, the latter reconstructs through a discommensuration which then moves towards the centre, 
where it arrives when the value of H is approximately that given by ([|). Next the discommensuration proceeds to 
broaden in discontinuous steps, with the AF-SF and SF-AF interfaces on the left and right side of the core moving 
outwards, until they arrive at the two surfaces when H is approximately at the bulk AF:SF transition value. The 
numerical results reported in 0,^,0] for even values of L > 16 and D — 0.5 agree well with the phase diagrams in 
Figs, [l] and § The same broadening scenario for the discommensuration, at a much smaller value of D, is visible in 
Fig 3 of §. 
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FIG. 1. Phase diagram in the H, D plane. Additional details are in Fig. 




FIG. 2. Phase diagram: details near the points P and Q. 
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